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Abstract 

The objective of this note is to present the results from the two papers [25] and [16]. We study the Navier- 
Stokes equation on the two-dimensional torus when forced by a finite dimensional white Gaussian noise. We give 
conditions under which both the law of the solution at any time t > 0, projected on a finite dimensional subspace, 
has a smooth density with respect to Lebesgue measure and the solution itself is ergodic. In particular, our results 
hold for specific choices of four dimensional white Gaussian noise. Under additional assumptions, we show that 
the preceding density is everywhere strictly positive. 

Resume 

Le but de cette Note est d'annoncer les resultats contenus dans les articles [25] et [16]. Nous etudions I'equation de 
Navier-Stokes sur le tore bidimensionnel forcee par un bruit blanc gaussien de dimension finie. Nous donnons des 
conditions sous lesquelles d'une part la loi de la solution a tout instant t > 0, projetee sur un espace de dimension 
finie, a une densite reguliere par rapport a la mesure de Lebesgue (qui, sous des hypotheses supplementaires, est 
strictement positive partout), et d'autre part la solution de la meme equation est un processus ergodique. En 
particulier ces resultats sont vrais dans certains cas de bruit blanc gaussien de dimension quatre. 



1. Introduction 



This note reports on recent progress made in [25,16] on the study of the two dimensional Navier-Stokes 
equation driven by an additive stochastic forcing. Recall that the Navier-Stokes equation describes the 
time evolution of an incompressible fluid. In vorticity form, it is given by 
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where x = {xi,X2) & T^, the two-dimensional torus [0, 27r] x [0, 27r], > is the viscosity constant, 
is a white-in-time stochastic forcing to be specified below, and 

B{w,w){x) = ^{Kw)i{x)^{x) , 



= 1 



where K. is the Biot-Savart integral operator which will be defined next. First, we define a convenient 
basis in which we will perform all explicit calculations. Setting = {(ji, j2) & : j2 > 0} U {(ji, ^2) G 

: ji > 0,j2 = 0}, Z2 = -Z% and Zl = Z\uZ^_, we define a real Fourier basis for functions on 
with zero spatial mean by 



ek{x) 



'sin(fc -x) e Z| 
cos(A; - a;) k & . 

Write w{t, x) = J2kei,^ otk{t)ek(x) for the expansion of the solution in this basis. With this notation, in 
the two-dimensional periodic setting, 

where = (— fc2,fci). See for example [20] for more details on the deterministic vorticity formulation 
in a periodic domain. We use the vorticity formulation for simplicity, but all of our results can easily 

be translated into statements about the velocity formulation of the problem. We solve (1) on the space 
= {/ = J2kez,2 "-kSk ■■ J2 lofcl^ < 00}. For / = J2kez,2 "-k^k, we define the norms = \ak\^ and 
\\f\\l = Y.\k?\ak?. 

The emphasis of this note will be on forcing which directly excites only a few degrees of freedom. Such 
forcing is both of primary modeling interest and is technically the most difficult. Specifically we consider 
forcing of the form 

W{t,x) = <^kWk{t)ek{x) . (3) 
fee 2, 

Here -Z* is a finite subset of Zq, cjk > 0, and {Wk ■ k G Z*} is a collection of mutually independent 
standard scalar Brownian Motions on a probability space (O, J^, P). 

The note describes two sets of results contained in two papers by two different subsets of the authors. In 
the first, Mattingly and Pardoux [25] give conditions ensuring that any projection of the time t transition 
probability of the solution of (1) onto a finite dimensional subspace has a C°° density with respect to 
Lebesgue measure. The result is based on the Malliavin calculus. Under additional conditions, this density 
is shown to be everywhere positive. The techniques developed are quite general and we expect they can 
be applied to many nonlinear, stochastic partial differential equations with additive noise. These results 
provide a first step towards a truly infinite dimensional version of Hormanders celebrated "sum of squares" 
theorem [17]. 

In the second paper, Hairer and Mattingly [16] give necessary and sufficient conditions for the main 
results and estimates of [25] to hold. They then use these tools to build a theory which, when applied 
to (1), proves that it has a unique invariant measure under extremely general and essentially sharp 
assumptions. In addition to the tools from [25], they introduce new concept and tool which together 
provide an abstract framework in which the ergodicity of (1) is proven. The concept is a generalization 
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of the strong Feller property for a Markov process which, for reasons that will be made clear below, is 
called the asymptotic strong Feller property. The main feature of this property is that a diffusion which 
is irreducible and asymptotically strong Feller can have at most one invariant measure. It thus yields a 
natural generalization of Doob's theorem. The tool is an approximate integration by parts formula, in the 
sense of Malliavin calculus, which is used to prove that the system enjoys the asymptotic strong Feller 
property. To the best of the authors knowledge, this paper is the first to prove ergodicity of a nonlinear 
stochastic partial differential equation (SPDE) under assumptions comparable to those assumed when 
studying finite dimensional stochastic differential equations. 

The ergodic theory of infinite dimensional stochastic systems, and SPDEs specifically, has been a topic 
of intense study over the last two decades. Until recently, the forcing was always assumed to be elliptic and 
spatially rough. In our context this translates to 2* = Zq and \ak\ ^ |^| " for some positive a. Flandoli 
and Maslowski [13] first proved ergodic results for (1) under such assumptions. This line of inquiry 
was extended and simplified in [11,14]. They represent a larger body of literature which characterizes 
the extent to which classical ideas developed for finite dimensional Markov processes apply to infinite 
dimensional processes. Principally they use tools from infinite dimensional stochastic analysis to prove 
that the processes are strong Feller in an appropriate topology and then deduce ergodicity. 

Next three groups of authors in [18,3,8], contemporaneously greatly expanded the cases known to be 
ergodic. They use the Foias-Prodi type reduction, first adapted to the stochastic setting in [21] and the 
pathwise contraction of the high spatial frequencies already used in [22] to prove ergodicity of (1) at 
sufficiently high viscosity. All of the results hinged on the observation that if all of the unstable directions 
are stochastically perturbed, then the system could be shown to be ergodic. A general overview of these 
ideas with simple examples can be found in [24]. These ideas have been continued in a number of papers. 
See for instance [7,4,23,15,19,24]. 

Unfortunately, the best current estimates on the number of unstable directions in (1) grow inversely 
with the viscosity i^. Hence the physically important limit of — > while a fixed, finite scale is forced 
were previously outside the scope of the theory. However there existed strong indications that ergodicity 
held in this case. Specifically in [9] it was shown that the generator of the diffusion associated to finite 
dimensional Galerkin approximations of (1) was hypoelliptic in the sense of Hormander when only a few 
directions were forced. This hypoellipticity is the crucial ingredient in the proof of ergodicity from [9]. 

The "correct" ergodic theorem needs to incorporate in its statement information on how the randomness 
spreads from the few forced directions to all of the unstable directions. This understanding when combined 
with what had been learned in [21,22,18,3,8] should yield unique ergodicity. This is the program executed 
in the papers discussed in this note. 



2. The Geometry of the Forcing and Cascade of Randomness 

The geometry of the forcing is encoded in the structure of from (3). As observed in [9], its structure 
gives information about how the randomness is spread throughout phase space by the nonlinearity. 

Define Zq to be the symmetric, and hence translationally stationary part of the forcing set 2* , given 
hy Zq = Zt, (1 {—Z^). Then define the collection 



Zr, = {e + jezl : jeZo, £eZn-i 



with i^-j^O, } 



and lastly, 



oo 
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2oo captures the directions to which the randomness has spread. This can be understood in the following 
way. Denote by dk the partial derivative into the direction efe of the phase space and define (on a formal 
level) the first order differential operator X by 

X= J2{B{w,w)k-iy\k\^)dk . 

Then the generator of the Markov process associated to (1) is formally given by 

C = X + ^Y.akdl. 

Note that B{w, w)k = J2e j '^k,j,eweWj, where Ck,j,e 7^ if and only if fc S {j ± £, —j ± £} and t'^ • j ^ 0, 
\j\ ^ \£\. Therefore, all differential operators of the type dk with k e Z^a can be obtained as an iterated 
Lie bracket of finite length involving Xo and with £ G . 

Since we want to ensure that all of the unstable directions are stochastically agitated, we seek conditions 
where Zoo — ^o- following essentially sharp characterization of this situation is given in [16]. 

Proposition 2.1 One has Z^o = if and only if both : 

(i) Integer linear combinations of elements of Zq generate Zq. 

(ii) There exist at least two elements in Zq with unequal euclidean norm. 

This characterization is sharp in the sense that if Z* = — and one of the above two conditions fails, 
then there exists a non-trivial subspace of which is left invariant under the dynamics of (1). Also 
notice that if 

Zo = {(0,1), (0,-1), (1,1), (-1,-1)} 

then Proposition 2.1 implies that Zoc = Zq. Hence forcing four well chosen modes is sufficient to have the 
randomness move through the entire system. Of course one can also force a small number of modes center 
elsewhere than at the origin and obtain the same effect. The next two sections discuss the implications 
of Zoo = Zg. 



3. Malliavin Calculus and Densities 

We define 

Soo = Span(efc : k G Zoo U Z*) . (4) 
One of the main results of [25] is the following : 

Theorem 3.1 For any i > and any finite dimensional subspace S of Soo, the law of the orthogonal 
projection I\.w{t, •) of w{t, •) onto S is absolutely continuous with respect to the Lebesgue measure on S 
and has a C°° density. 

In [10], Eckmann and Hairer used Malliavin calculus to prove a version of Hormandcr's "sum of squares" 
theorem for a particular SPDE and deduce ergodicity. However, all of the techniques of that paper required 
that the forcing excite all but a finite number of directions and that the forcing be spatially rough as in 
[13,6]. The proof of Theorem 3.1 builds on ideas introduced into Malliavin calculus by Ocone in [26]. The 
central idea is an alternative representation of the Malliavin matrix of (1) using the time reversed adjoint 
of the linearization of (1). Ocone used this representation when the SPDE was linear in the initial data 
and the forcing. When the noise is additive, [25] extends that idea to the nonlinear case. 

Let Js,t^ be the solution of linearization of (1) at time t with initial condition ^ at time s, s < t. Let 
J*(^ denote the solution to the L^-adjoint of the linearizion at time s, s < t, with terminal condition ^ 
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at time t. Since the equation is time reversed, the adjoint is well posed. With this notation, the so-called 
"Malliavin covariance matrix" Ait can be represented by 




where (/> G L^. The second of these representations is the one used in [25]. Because of the time reversal, the 
representation is not adapted to the filtration generated by W and new estimates concerning anticipating 
stochastic processes are required to obtain the needed estimates. Essentially one needs to show that the 
Malliavin matrix is non-degoncratc on the subspacc S and that the moments of the reciprocal of the norm 
of the Malliavin matrix on this subspace are finite. This is accomplished through the following estimate 
which also gives information about the separation of the randomness on large and small scales. 

Theorem 3.2 Let U be the orthogonal projection ofL? onto a finite dimensional subspace of Soo- For 
any t > 0, -q > 0, p > 1, M > and K G (0, 1) there exist two constants c = c{y,r],p, \Z^\,t,K,M,Yl) 
and eo = eo{u,K, |Z*|,t,M,n) such that for all e G (0,eo]; 



With additional assumptions on the controllability of (1) conditions are also given ensuring the strict 

positivity of the density. This extends results of Ben Arous and Lcandre [2] and Aida, Kusuoka and 
Stroock [1] to this setting. We refer the reader to [25] for the exact conditions and the details. 



4. Unique Ergodicity 

Recall that an invariant measure for (1) is a probability measure /x* on such that Ptfi* — /i^, where 
Pj* is the semigroup on measures dual to the Markov transition semigroup Pt defined by {Pt(t>){w) = 
E,yj(j){wt) with (j) G Cb{h'^). While the existence of an invariant measure for (1) can be proved by "soft" 
techniques using the regularizing and dissipativity properties of the flow [5,12], showing its uniqueness 
is a more challenging problem that requires a detailed analysis of the nonlinearity. The importance of 
showing the uniqueness of /x* is illustrated by the fact that it implies that 



for all bounded continuous functions and all initial conditions wq G . It thus gives some mathematical 
ground to the ergodic assumption usually made in the physics literature when discussing the qualitative 
behavior of (1). The following theorem is the main result of [16]. 

Theorem 4.1 If = Zg, then, (1) has a unique invariant measure in L^. 

When combined with Proposition 2.1, this theorem gives easy to verify conditions guaranteeing a unique 
invariant measure. 

The concept of a strong Feller Markov process appears to be less useful in infinite dimensions than in 

finite dimensions. In particular if Pi is strong Feller, then the measures Pt{u, •) and Pj(t', ■) arc equivalent 
for all initial conditions u,v £ h"^ . It is easy to construct an ergodic SPDE which does not satisfy this 
property. 




where S{M,K,Il) = {(P G : ||<^|| = 1, ||</.||i < M, \\Il<f>\\ > K}. 




(5) 
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Recall the following standard sufficient criteria for Pf to be strong Feller : there exists a locally bounded 
function C{w,t) such that 

\V{Pt4>){w)\<C{w,t)\\4>\\oo 

for all Frechet differentiable functions : M. While wc will not give the exact definition of the 

asymptotic strong Feller property here, the following similar condition implies that the process is asymp- 
totically strong Feller : there exists a locally bounded C{w), a non-decreasing sequence of times tn, and 
a strictly decreasing sequence e„ with e„ ^ so that 

|V(Pt„</.)H| <CH||.^||oo+e„||V.^||oo (6) 

for all Frcchct differentiable functions </) : ^ R and all n > 1. In applications one typically has i„ ^ oo. 
Hence, the process behaves as if it acquired the strong Feller property at time infinity, which justifies the 
term asymptotic strong Feller. 

First observe that (V^(Pt0)(w), ^) = Et„(V(/))(wt) Jo.f^. Next wc seek a direction v in the Cameron- 
Martin space so that if denotes the Malliavin derivative in the direction v then Jo,fC = T>^wt. In finite 
dimensions, we can often do this exactly ; however, in infinite dimensions we only know how to achieve 
this up to some error. Setting pt = Jo,tC ~ T^^wt, we have the approximate integration by parts formula. 

E^(V(^)K) Jo,t^ = + E^(V<^)K)Pt 

Jo 

From this equality one can quickly deduce (6), provided E| VsdWs\ < oo and E|pt| ^ as f ^ oo. 
In [16], a Vt is chosen so that these conditions hold. The analysis is complicated by the fact that the vt 
constructed there is not adapted to the Brownian filtration. This complication seems unavoidable. Hence, 
the stochastic integral is a Skorohod integral and all of the calculations are made more complicated. 

The ideas developed here can also be used to prove exponential mixing using the ideas from [23,15]. 
These results will be presented elsewhere. 
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